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We study the problem of whether the active gravitational mass of an isolated sys- 
tem is equal to the total energy in the tetrad theory of gravitation. The superpotential 
is derived using the gravitational Lagrangian which is invariant under parity opera- 
tion, and applied to an exact spherically symmetric solution. Its associated energy is 
found equal to the gravitational mass. The field equation in vacuum is also solved 
at far distances under the assumption of spherical symmetry. Using the most general 
expression for parallel vector fields with spherical symmetry, we find that the equality 
between the gravitational mass and the energy is always true if the parameters of the 
theory a±, a<i and a 3 satisfy the condition, {a\ -\-a-i){a\ — ^ 0. In the two special 

cases where either (a± + a?) or (a± — is vanishing, however, this equality is not 

satisfied for the solutions when some components of the parallel vector fields tend to 
zero as 1/ 'y/r for large r. 
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§1. Introduction 



It was shown by Miller 1 ) that a tetrad description of gravitational field allows a more 
satisfactory treatment of the energy-momentum complex than general relativity (GR). 
The Lagrangian formulation of the theory was given by Pellegrini and Plebanski. 2 -* 
Hayashi and Nakano 3 -* independently formulated the tetrad theory of gravitation as a 
gauge theory of the spacetime translation group. In the earlier attempts, admissible La- 
grangians were limited by the assumption that the field equation has the Schwarzschild 
solution. Miller 4 -* later suggested to abandon this assumption and to look for a wider 
class of Lagrangians. Meyer 5 -* formulated the tetrad theory as a special case of Poincare 
gauge theory. 6 -*' 7 -* Saez 8 -* generalized the theory into a scalar-tetrad theory. 

Hayashi and Shirafuji 9 -* studied the geometrical and observational basis of the tetrad 
theory of gravitation. Geometrically the tetrad fields are identified with the parallel 
vector fields defined by the underlying absolute parallelism. Incidentally they gave the 
name, new general relativity (NGR), to the theory of gravitation based on absolute 
parallelism, since Einstein 10 - 1 was the first to introduce the notion of absolute parallelism 
in physics. They assumed the Lagrangian with three unknown parameters, denoted 
by ai, a 2 and a 3 . In order to reproduce the correct Newtonian limit, the first two 
parameters should satisfy a condition called the Newtonian approximation condition, 
which allows us to express these two parameters in terms of an unknown dimensionless 
parameter e. An exact, vacuum solution of the gravitational field equation was found for 
static, spherically symmetric case, where the parallel vector fields take a diagonal form. 
This solution describes the spacetime around a mass point located at the origin, and will 
be referred to as the "exact spherically symmetric solution" hereafter. Then comparison 
with solar-system experiments showed that the value of e should be very small. By 
contrast only an upper bound has been estimated for the remaining parameter a^. 9 ^' 11 ^ 
The singularity problem of the exact solution has been studied. 12) ' 13) If (ai + a 2 )=0, 
then the theory reduces to the one studied by Hayashi and Nakano 3 - 1 and Miller. 4 -* 

The equation of motion for a test particle was discussed based on the response 
equation of the energy-momentum tensor of matter fields. 9 -* In particular, when the 
intrinsic spin of the fundamental particles constituting the test particle can be ignored, 
the response equation reduces to the covariant conservation law of GR, and hence the 
world line of the test particle is a geodesies of the metric defined by the parallel vector 
fields. Accordingly, for spinless test particles the passive gravitational mass is equal to 
the inertial mass. 

In GR the equality between passive gravitational mass and inertial mass is pre- 
dicted to be valid also for massive bodies like planets which contains appreciable frac- 
tion (about 5 x 1CT 10 for the Earth) of gravitational self-energy. This equality means 
that gravity pulls on gravitational binding energy of a massive body just as it does on 
other forms of mass-energy, and is referred to as the strong equivalence principle. 14 -* 
Nordtvedt 15 -* developed the parametrized post-Newtonian (PPN) formalism and calcu- 
lated the diviation from the unity of the ratio of gravitational mass to inertial mass 
in terms of the PPN parameters. Analysis of the lunar laser ranging (LLR) over past 
24 years has confirmed that the Earth and the Moon accelerate equally to the sun 
within fractional difference less than 5 x 10~ 13 , which yields the Nordtvedt parame- 
ter i] = —0.0010 ± 0.0010. 16 )' 17 )' 18 ) This parameter i] (= in GR) can be expressed 
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as r] = 4/5 — 7 — 3, in fully conservative theory which possesses a full complement 
of the post-Newtonian conservation laws: energy, momentum, angular momentum 
and center-of-mass motion. 14 ) Here (3 and 7 are the so-called Eddington- Robertson 
parameters. 19 ) ,2 °) 

Thus, if NGR is a fully conservative theory in the above sense, then the Nordtvedt 
parameter 77 vanishes, since the Eddington- Robertson parameters are expressed as j3 = 
1 — e/2 and 7 = 1 — 2e in terms of the dimensionless parameter e mentioned above. It 
can be shown that energy, momentum and angular momentum are conserved in NGR. 
We must be careful not to conclude that rj = 0, however, since NGR is not a metric 
theory: Accordingly we are now trying to develop its PPN formalism. 

It is the purpose of this paper to study a different aspect of the equivalence prin- 
ciple, namely the problem of whether or not the active gravitational mass (or simply 
the gravitational mass) of an isolated system is equal to its inertial mass, i.e. the total 
energy divided by the square of the velocity of light As is well known, this problem is 
settled affirmatively in GR. 21 -* In the case of (a\ + 02) = Mikhail et al. 22 -* calculated 
the energy of two spherically symmetric solutions, and found that the energy in one 
of the two solutions does not coincide with the gravitational mass. Shirafuji et al. 23 -* 
extended the calculation to all the stationary asymptotically flat solutions with spher- 
ical symmetry, dividing them into two classes: The one in which the components, (& a o) 
and (b^°' a ), of the parallel vector fields (b k At ) tend to zero faster than 1/y/r for large r 
and the other in which those components go to zero like 1 / y/r .f\ It was found that the 
equality of the gravitational and inertial masses holds true only in the first class. 

In this paper we study the problem for stationary, spherically symmetric systems 
without making any assumption on the parameters, ai, a 2 and a 3 other than the 
Newtonian approximation condition. 

We organize this paper as follows. In §2 we give a brief review of the NGR. In §3 we 
derive the superpotential of the total energy-momentum complex, and apply it to the 
exact spherically symmetric solution. In §4 we study stationary spherically symmetric 
solutions at far distances in the linear approximation. In §5 we calculate the energy in 
the special case of (a\ — Aas/9) = 0, taking all the leading terms into account beyond 
the linear approximation. The final section is devoted to conclusion and discussion. 



§2. A brief review of the NGR 



We assume spacetime to admit absolute parallelism, i.e. to have a quadruplet of 
linearly independent parallel vector fields (b k At ) satisfying 

djK = b\,u - r V"a = (i) 

* Since we use the unit c = 1, we need not draw a distinction between the inertial mass and the 
total energy 

tin this paper Latin indices (i,j,k, ...) represent the vector number, which runs from (0) to (3), 
while Greek indices (/x, v, p, ...) represent the world-vector components running from to 3. The 
spatial part of Latin indices are denoted by (a, b, c, ...), while that of Greek indices by (a, (3, 7, ...). 
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with & fc A t,i/ = d v b k fl . Solving (1), we get the nonsymmetric connection 

r A M i/ = bk X b k ^j,, (2) 

which defines the torsion tensor as 

T ^ = V ^ v r up = 6^ (b ^ v b i/,n) ■ (3) 

The curvature tensor defined by T x ^ u is identically vanishing, however. The metric 
tensor is given by the parallel vector fields as 

9»» = bk^b k U} (4) 

where we raise or lower Latin indices by the Minkowski metric Tjij—r]^ =diag(-l,+l,+l,+l). 

Assuming the invariance under 

a) the group of general coordinate transformations, 

b) the group of global Lorentz transformations, and 

c) the parity operation, 

we write the gravitational Lagrangian density in the form [] 



r - V-9 

*-G — 



ai + a 2 (v^) + a 3 (a^a^)} , (5) 



where ai, a 2 and a 3 are dimensionless parameters of the theory, f\ and 

Va = \ ( T ^x + T UfiX ) + ^ {gx^u + gxuv^) - ^g^vx, (6) 

v, = T\„ (7) 

(III "^^uupaT ^ (8) 

o 

with e^ upa being the completely antisymmetric tensor normalized as e i23 = \J—Q- By 
applying variational principle to the above Lagrangian, we get the field equation: 

jn* = kT i»> ( 9 ) 

with 

1^ = 2k[D x F^ uX + v x F^ x + - -g^Lo], (10) 



* Throughout this paper we use the relativistic units, c = G = 1. The Einstein constant k is then 
equal to 87r. We will denote the symmetric part by the parenthesis ( ) and the antisymmetric part by 
the square bracket [ ] . 

' The dimensionless parameters na,i of Ref.9) are here denoted by ai for convenience. 



4 



where 



W v\ = 1 h kn dL G 
~ 2 dbK,x 



L G 

rpflV 



1, 



rppafj, p v _rpvpa ff up 



1 SC M 
/= g 5b k u 



(11) 
(12) 

(13) 
(14) 



Here Cm denotes the Lagrangian density of material fields and T pv is the material 
energy-momentum tensor which is nonsymmetric in general. In order to reproduce the 
correct Newtonian limit, we require the parameters d\ and a 2 to satisfy the condition 



ai + 4a 2 + 9axa 2 = 



called the Newtonian approximation condition, 9 -* which can be solved to give 



(15) 



1 



ai = 



3(1 -eV 



a 2 = 



3(1 -4e) 



(16) 



with e being a dimensionless parameter. The comparison with solar-system experiments 
shows that e should be given by 9 ) 



e = -0.004 ±0.004, 



(17) 



which we assume throughout this paper. 

It is well known that the conservation law in GR is given by 



wv = o, 



(18) 



where Tq^ v is the symmetric material energy-momentum tensor of GR and the semi- 
colon denotes covariant derivative with respect to the Christoffel symbol. This law 
does not follow from (9), however. Instead, we can derive the response equation 



T^, u = K^T [uX] , 



(19) 



where K uXfM is the contortion tensor given by 



p-uXfl ^_ ( rpl/\[l p\vfl rptiv\\ J^XvjX 

~ 9 V ) ~ 



(20) 
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The antisymmetric part is due to the contribution from the intrinsic spin of 

fundamental spin-1/2 particles. For macroscopic test particles for which the effects 
due to intrinsic spin can be ignored, their energy-momentum tensor can be supposed 
to be symmetric and satisfy (18). The equation of motion for such macroscopic test 
particles is then the geodesic equation of the metric. 

In static, spherically symmetric spacetime the parallel vector fields take a diagonal 
form, and the field equation (9) can be exactly solved in vacuum to give 



' m " 




pr 






q/2 


(»-) 




V qr) 










\ 



\ (2-p)/2 / 

m_] 1 + m 



\ (2+?)/2 



pr ) 



qr ) 



(21) 



a J 



where m is a constant of integration, and the constants p and q are given by 



p 



■[(1 -5e + 4e 2 )2 - 2e], q 



(l-5e) LV ' J ' ^ (l-5e) 

In spherical polar coordinates (21) gives the line-element 



[(1 -5e + 4e 2 )2 +2e] 



(22) 



ds 2 = 



1 + qr) 



)p 
^dt 2 +(l-™ 
m\ V pr 



2-p 



qr ) 



2+q 



dr 2 + r 2 (dd 2 + sin 2 dd<f 



(23) 



^From the asymptotic behavior of the component goo of the metric tensor, the constant 
m can be identified with the gravitational mass of the central gravitating system. 
It is clear that if the parameter e vanishes, the line-element (23) coincides with the 
Schwarzschild solution written in the isotropic coordinates. 



§3. Superpotential of the NGR and calculation of the energy 



The following identity can be derived from the invariance of Cq under general coor- 
dinate transformations ■}^ 24 ^ 



When (6%) satisfies the gravitational field equation (9), this implies 

(T/ + */) = dx (2v^/ A ) , 
6 



(24) 



(25) 



where is the canonical energy-momentum complex of the gravitational field 



uk , r v r 



(26) 



The total energy-momentum complex is then defined by 



M, 



(27) 



with IA^ X being the superpotential 



W/ A = 2^F M 



(2* 



Since the tensor F M I/A is antisymmetric with respect to v and A, the M.^ of (27) satisfies 
ordinary conservation law, 

(29) 



d v M; = 0. 



The total energy is now given by 



E — — [ M °d 3 x = - lim / 



r=constant 



U 0a n a dS, 



(30) 



where n a is the outward unit 3-vector normal to the surface element dS. 

Let us calculate the superpotential by writing the Lagrangian (5) in the form 



C 



'-9 



G 



(31) 



where L« = t x ^t x ^, £ (2) = v^Vn and = a^a^. Writing (11) in the form 



uX 



uX 



a i 



uX 



(32) 



we get 



i?(2) 

i?(3) 



uX 1 



vX 



vX 



2T/ A + T x \ - T uX ^ - (^V - V) 

(<^V - V) > 

1 

~9 



t/ x -t x % + t 



-\uX 



(33) 
(34) 
(35) 



*The Lagrangian used by Miller 4 -* is different from that used by Hayashi and Shirafuji 9 -' by a factor 
(-2). Accordingly, the definition (28) is different from that of Miller by a factor (-2). 
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where F^^, F ( - 2 ^ uX and F^ ^ X correspond to and L^ 3 \ respectively. So 

with the help of (28) the superpotential of the NGR can be written as 



2v^ 



K 



a 3\r rV \ . f a l . a 3\ (rpXv 



9 



rpv X \ 
1 N 



a i 



a 2 ) (^V-^V) 



(36) 



As an example, let us apply (36) to the exact solution (21). The appropriate 
components U 0a are given by 



m n 
2npq r 2 



TTX 

(p-q- 8)(ai - 2a 2 )— + 4(p - q)(ai - 2a 2 ) 

r 



+Apq(2 ai - a 2 ) + 3(p - q) 



aim 



(37) 



with p and g given by (22). Using (37) in (30), we get 



E = 



8irm 
npq 



[{p ~ 0) (°i - 2fl 2) + pg(2ai - a 2 )] 



= m [(2ai — a 2 ) — 2e(ai — 2a 2 )] = m, 



(38) 



where the relation {p — q)/pq = — 2e, which follows from (22), is used in the second 
equality and (16) is employed in the last one. This means that the total energy of 
the exact solution is just the same as the gravitational mass of the central gravitating 
system. 



§4. The spherically symmetric parallel vector fields and its 
energy in the linear approximation 

The solution (21) is the only exact spherically symmetric solution in vacuum that 
we know at present, and it is not clear to us whether there exist any other spheri- 
cally symmetric solutions in vacuum. In view of this let us consider a wider class of 
spherically symmetric solutions at far distances from the source. Consider an isolated, 
gravitating system with spherical symmetry, and restrict attention to the weak field 
far from the source. The most general form of the parallel vector fields is given in the 
Cartesian coordinates by 25 -* 



C{r) G{r)n a 
H(r)n a 5 a a D(r) + E(r)n a n a + F(r)e (lQ ^n /3 



(39) 
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where the two real functions C(r) and D(r) are supposed to approach 1 at infinity, 
while the remaining four real functions, E(r), F(r), G(r) and H(r), must tend to zero 
there. Here we define the radial unit vector n a and n a by 



n 



.r" 
r 



5 a a n\ 



(40) 



without making distinction between upper and lower indices. Using the freedom to 
redefine the radial coordinate r, we can eliminate the function E(r) from the compo- 
nents (b a a ). Accordingly we can put E(r) = without loss of generality. The metric 
tensor is then written as 



g 00 = -(C 2 - H 2 ), 

9oa = {-CG + DH}n a , 

g aP = (D 2 + F 2 )5 aP - (F 2 + G 2 )n a n p . 



(41) 



According to (30), the total energy of an isolated system can be calculated if the 
superpotential is known up to order 0(l/r 2 ). It is then enough to know the parallel 
vector fields up to order 0(l/r). So let us restrict our attention to the weak field far 
from the source, and analyze the field equation in vacuum up to order 0(l/r 3 ). This 
has been performed in GR, 21 -* showing quite generally that the gravitational mass is 
equal to the total energy for any stationary isolated system. In the NGR, however, the 
analysis of the field equation in vacuum has not yet been done for the general case of 
(ai + a 2 ) 7^ even in the weak field approximation. This is due to the fact that the 
gravitational field equation of the NGR is more complicated than the Einstein equation 
of GR. 

Let us suppose that the leading term of the five unknown functions is given by some 
power of 1/r, and that (b k ) can be represented as f\ 



v 



-n 



\ 



(42) 



Here the powers, s, t, u, v and w, are positive unknown constants at the beginning 
of the calculation, and their value will be determined by the field equation of NGR. 
Here the constant coefficients, b, j, h, d and /, are also unknown, but they can be 
assumed to be nonvanishing without loss of generality since the powers of r are left 
unknown. Use (42) in (41) gives the asymptotic behavior of the metric tensor which 
involves linear and quadratic terms of the unknown constants, b, j, h, d and /: The 
linear terms are dominant if the powers satisfy the following inequality 



s < 2u, 



v < 2t, 



v < 2w, 



(43) 



The constant b should not be confused with det(6 fc M ), which we denote by 
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which we call the condition of the linear approximation. 

Now apply the vacuum field equation (9) to (42), assuming the condition (43). 
Keeping only the leading terms, which are shown to be linear in the five unknown 
constants due to (43), we get the nonvanishing components of I^ u ] 

-2 (s(s- l)(ai + a 2 )b v(v - l)( ai - 2a 2 )d\ 
^o-^j ~, ~v J> ^ 

I a0 = 2(u + l)(u-2) {ai+ u+ a 2 2)h n a: (45) 



'a/3 



n a n/3 \ s(s + 2)(ai — 2a 2 )b v(v + 2)(ai + 4a 2 )d 



/ s 2 (ai — 2a 2 )b f 2 (ai + 4a 2 )d 

2 i y S ryl) 



r 



+(w + l)(w - 2)^^ e aPl n\ (46) 

Here it is important to notice that I a o and I[ a p] are dominated by the linear term 
irrespectively of the condition of the linear approximation. Also using (42) in (36), we 
get the leading terms for the appropriate components of the superpotential; 



9n a 



s(ai + a 2 )b v(ai — 2a 2 )d 



(47) 



By virtue of (44) and (46) we can show that the equations, Joo = and I^ a p) = 0, 
give the following results: 

1) if s < v then 6 = 0, because (ai — 2a 2 ) ^ 0, 

2) if s > v then d = 0, because (ai + 4a 2 ) 7^ 0, and 

3) if s — v 7^ 1 then b — d — 0, because 

/ (ai + a 2 ) -(a 1 -2a 2 )\ 
det ^ 0. (48) 

\(a l -2a 2 ) -(ai + Aa 2 )J 

All the above three cases contradict our assumption that b 7^ and d 7^ 0. Therefore, 
we find that the powers s and v should be given by 

s = v = l. (49) 

Then the equation 7 00 = is automatically satisfied and the remaining one, I{ap) — 0, 
gives 

d = (2e- 1)6. (50) 
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We notice that the asymptotic behavior of the diagonal exact solution (21) satisfy (49) 
and (50) with b = —m. Also this result of (49) is physically very satisfactory in view 
of the superpotential (47), because we can then get a finite value of the energy. 

Next the (aO)-component of the field equation in vacuum, I a0 = 0, implies that 

(u-2)( ai + a 2 ) =0. (51) 

When (ai + a 2 ) ^ 0, this gives 

u = 2, (52) 

compatibly with the condition of the linear approximation (43). Using (49) and (51) 
in (41), we see that the component g 00 of the metric tensor behaves asymptotically like 

9oo = - U + -), (53) 



which indicates that the constant c is related to the gravitational mass m of the isolated 
system by 

b = -to, (54) 

since as shown in §2 the world line of a spinless test particle is a geodesies of the metric. 
So (50) can now be written as 

d = (1 - 2e)m. (55) 

In the special case of {a\ + a 2 ) = 0, on the other hand, the power u is not constrained 
by the field equation I a0 = 0. However, the exact form of all the spherically sym- 
metric solutions was found in this special case, and the energy of those solutions was 
calculated, 23 -* as will be explained at the end of this section. 

Finally from the skew part, I[ a ^ = 0, we get 

(w - 2) ( fll - ^) = 0. (56) 

When (ai — 4as/9) ^ 0, it follows from the condition (56) that 

w = 2, (57) 

satisfying the condition of the linear approximation (43). In the special case of 

(ai — 4a 3 /9) = 0, on the other hand, the power w is not restricted by the field equation. 

We shall discuss this case separately in the next section. 

Collecting the above arguments together, we see that when the parameters satisfy 
(ai + a 2 ){a\ — 4a3/9) 7^ 0, the field equation = can be solved in the linear 
approximation to give the following asymptotic form of (b k fJi ): 



/a- 



m n a 
7) 



\ hy 2 
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(b% (exact)) + r* + - , (58) 



where [b ^ (exact) J denotes the exact solution (21) in the asymptotic form. Substitut- 
ing (58) into (36) gives 



Ko° a = r [(ai + a 2 ) + (1 - 2e)( ai - 2a 2 )] = (59) 



where we have used (16) in the last equation. From (30) we then get 

E = m. (60) 

Thus, the gravitational mass is equal to the total energy for an isolated spherically 
symmetric system. 

For completeness let us recapitulate the known results 23 ) about the exact, spheri- 
cally symmetric solutions and their total energy. Referring to the general expression 
for the parallel vector fields given by (39), the solutions are divided into two classes: (1) 
the solution with F(r) = and (2) the solution with F(r) ^ 0. The general solution 
of the class (1) involves an arbitrary function of r, and therefore there exist solutions 
whose components (b a ) asymptotically behave like l/r u for any positive value of u. 
The calculated energy was shown to coincide with the gravitational mass only when 
the components (6 a o) go to zero faster than l/yjr. When (6 a o) ~ hn a / 'y/r for large 
r, the quadratic terms of h must be taken into account in the superpotential, and the 
energy does not coincide with the gravitational mass, being given by 

h 2 

E = m+—. (61) 

If (6 a o) behave like l/r u with < u < 1/2, the calculated energy is infinite, implying 
that the solutions with such an asymptotic behavior is physically unacceptable. Next 
let us turn to the general solution of the class (2), which was shown to involve an 
arbitrary constant parameter, and to have the asymptotic behavior of (42) with h = 
and w = 2: The calculated energy was shown to coincide with the gravitational mass. 



§5. The energy in the special case of (a\ — 4<23/9) = 

Now let us turn to the special case of (ai — = 0, in which we must go beyond 

the linear approximation discussed in the preceding section since the field equation 
does not impose any restriction on the power w. We notice that when the parameters, 
ai, a 2 and a 3 satisfy the conditions, (ai + a 2 ) = and (ai — 4a 3 /9) = 0, together with 
the Newtonian approximation condition (15), the gravitational Lagrangian (5) reduces 
to the Einstein-Hilbert Lagrangian of GR written in terms of the tetrad field, thus 
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leading to inconsistency of the field equation (9): For example, the is symmetric, 
while T^ u of spinor fields is nonsymmetric. Accordingly we assume {a\ + a 2 ) 7^ in 
this special case. 

Now consider the parallel vector fields (b k ^) which asymptotically behave like (42) 
with u = 2. Since the power w cannot be fixed by the field equation, we shall discuss 
the three cases, w > 1/2, w — 1/2 and w < 1/2, separately. 

(i) When w > 1/2, it can be shown by the same argument as in the previous section 
that the powers s and v satisfy either s — v — 1, or s > 1 and v > 1. If s — v — 1, 
the linear approximation studied in the preceding section is still valid, and the parallel 
vector fields behave asymptotically like 



(b\) = (b\{exact)) + 



/0 \ 

f 



(62) 



and the energy coincides with the gravitational mass. 

If it happens that s > 1 and v > 1, on the other hand, the superpotential dies 
out at infinity faster than 1/r 2 , and hence the calculated energy will be vanishing. 
The gravitational mass is also vanishing since goo tends to 1 for large r faster than 
1/r. Therefore, such a solution, if it exists, is devoid of physical meaning, although we 
cannot exclude its existence by the present approximate treatment. 

(ii) When w = 1/2, the field equation I af3 = implies that the powers s and v must 
satisfy one of the three possible alternatives: (a) s — v — 1, (b) s — 1 and v > 1 or (c) 
s > 1 and v — 1. 

Let us start with the case (a) of s — v — 1. The left-hand side of (9) is given by 



'00 — J-Oa — J-aO 



0. 



lap = 



6(ai - 2a 2 ) - d{ai + 4a 2 ) + 2/ 2 (ai + a 2 ) 



'8 a p - 3n a np 



(63) 



up to order 0(l/r 3 ), where we have used (ai — 4as/9) = 0, and the constant c satisfies 
(54) since s — 1 and u — 2. The vacuum field equation then provides the condition 



d = 



a\ + 4a 2 



6(ai - 2a 2 ) + 2/ 2 (ai + a 2 ) = (1 - 2e)m + 2e/ 



(64) 



Similarly the superpotential is expressed by 

2n a 



0a 



1 



ai + a 2 )6 - (ai - 2a 2 )rf + ^(ai - 2a 2 )/ 2 
/tr z L' 2 

2mn Q r 1 - 2e / 2 " 
1 — 

Kr 2 1 — e 2m ' 



(65) 
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where we have used (16) as well as (64). With the help of (30) we then have 



In the case (b) of s = 1 and v > 1, the terms proportional to the constant d do not 
contribute to the leading term of I a p and U 0a . Accordingly, the calculated energy is 
written as 

which goes to infinity as e tends to zero (i.e., (a± + a 2 ) — * ). This means that when 
the parameter e is varied, the solution of the case (b), if it exists, will be singular at 
e = 0. 

The last case (c) of s > 1 and v — 1 can be examined in the similar manner. Since 
s > 1, the gravitational mass of the central body must be vanishing, and furthermore 
the energy takes a value 

E = -^^-j\ (68) 
2(1 -e) J ' V ; 

which is negative due to the experimental value of e given by (17). Thus, the case (c) 
should be discarded. 

Therefore, when the power w is equal to 1/2, the energy of the spherically symmetric 
body differs from the gravitational mass as measured from far distances. It must be 
noticed, however, that we have not yet actually found a solution with such asymptotic 
behavior with w = 1/2. The above result merely implies that the existence of such a 
solution is not excluded by the field equation which considers only up to order 0(l/r 3 ): 
Detailed analysis based on exact solutions is desirable. 

(iii) When w < 1/2, quadratic terms of / contribute to the superpotential with the 
asymptotic behavior like l/r 2w , and hence the energy integral will diverge, giving an 
infinite value. Thus, the solution with w < 1/2, if it exists, must be discarded. 



6. Conclusion and discussion 



We studied the problem whether or not the equality of the active gravitational 
mass and the inertial mass (i.e., the total energy) is ensured by the gravitational 
field equation in NGR, restricting ourselves to isolated spherically symmetric systems. 
Based on the identity following from the general coordinate invariance, we defined the 
gravitational energy-momentum complex and derived the explicit expression for the 
superpotential which allows us to calculate the total energy of isolated systems. We 
first applied it to the exact spherically symmetric solution of diagonal matrix form, 
and found that the equality of the active gravitational mass and the total energy is 
satisfied. 
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Since we do not know whether the spherically symmetric solution is unique, we then 
discussed the asymptotic method to calculate the total energy. Namely, we started from 
a quite general expression for the parallel vector fields with spherical symmetry, and 
solved the vacuum field equation at far distances from the source. The main results 
can be summarized as follows: 

(a) In the generic case of (ai + a 2 )(ai — 4a 3 /9) 7^ 0, the linear approximation can be 
applied to solve the field equation. All the components of the parallel vector fields 
other than (b^ a ) are asymptotically determined by the field equation, and coincide 
with those of the exact solution up to order 0(l/r). The components (b^ a ) do not 
contribute to the total energy, however, and the calculated value of the energy agrees 
with the gravitational mass. 

(b) In the special case of (a 1 + a 2 ) = the asymptotic behavior of the components (b a ) 
is not fixed at all by the field equation. We know, however, that for any asymptotic 
behavior of (6%), there exist exact solutions with such asymptotic behavior. 23 ) The 
total energy is finite only when (&%) ~ l/r u with u > 1/2, and furthermore if u = 1/2, 
the quadratic terms of (b a ) contribute to the total energy, violating the equality of the 
total energy and the gravitational mass. 

(c) In the special case of (ai — 4a 3 /9) = 0, the gravitational field equation does not 
impose any condition on the e aa/ gn^-term of (b a a ). We showed that if this term tends 
to zero at infinity like 1/y/r, the equality between the gravitational mass and the total 
energy is violated. It is yet to be studied, however, whether exact solutions with such 
asymptotic behavior really exist. 

Thus, in the special cases of (b) and (c) above, there are anomalous solutions which 
violate the equality of the energy and the gravitational mass. The characteristic feature 
of these solutions is that specific components of the parallel vector fields tend to zero as 
l/\fr for large r. The physical meaning of these anomalous solutions is not yet clear, 
however. Is it reasonable to rule out such anomalous solutions by demanding that all 
the components of the parallel vector fields should tend to flat-space value faster than 

As we noticed, if the parameters satisfy both the conditions, (01 + 02) = and [a\ — 
4a 3 /9) = 0, then the gravitational Lagrangian in the NGR coincides with the Einstein- 
Hilbert Lagrangian in GR expressed in terms of the tetrad field. The gravitational 
field equation of the NGR then becomes inconsistent: In fact, the L.H.S. of (9) will be 
symmetric but the R.H.S. is nonsymmetric. Therefore, we assumed in this paper that 
the combinations of the parameters, (01+02) and (a\ — 4a 3 /9), do not vanish at the same 
time. The special case of (ai + a 2 ) = has been studied in rather detail. 3 )+' 9 ): 22 )> 23 ) 
Another case of (ai — 4a 3 /9) = deserves more attention in this respect. 

We assumed the spherical symmetry to obtain the solution at large distances. In 
GR, however, it is well known that such assumption of spherical symmetry can be 
removed. In fact, Einstein equation ensures that the metric tensor at far distances is 
the same as the Schwarzschild metric up to order 0(1/ r) for any isolated stationary 
system. Is it possible in NGR to investigate the vacuum solution at far distances 
without assuming spherical symmetry? 

We have ignored the parity- violating term, i> /1 a M , throughout this paper. The pos- 
sibility that this term plays an important role was suggested by Muller-Hoissen and 
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Nitsch 26 ** in the case of (ai + a 2 ) = 0. It seems interesting to take into account this 
parity-violating term also in the case of {a\ — 4a 3 /9) = 0. 

We shall address ourselves to these problems in future work. 
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